Phase diagram of the XXZ ferrimagnetic spin-(l/2, 1) chain in the presence of 

transverse magnetic field 



A. Langari-'^, J. Abouie^-'^, M. Z. Asadzadeh^ and M. Rezai-*^ 
^ Department of Physics, Sharif University of Technology, Tehran 11155-9161, Iran 
^ Department of Physics, Shahrood University of Technology, Shahrood 36199-95161, Iran and 
School of Physics, Institute for Research in Fundamental Sciences (IPM), Tehran 19395-5531, Iran 

(Dated: January 14, 2013) 

We investigate the phase diagram of an anisotropic ferrimagnetic spin-(l/2, 1) in the presence of a 
non-commuting (transverse) magnetic field. We find a magnetization plateau for the isotropic case 
while there is no plateau for the anisotropic ferrimagnet. The magnetization plateau can appear only 
when the Hamiltonian has the U(l) symmetry in the presence of the magnetic field. The anisotropic 
model is driven by the magnetic field from the Neel phase for low fields to the spin-fiop phase for 
\ intermediate fields and then to the paramagnetic phase for high fields. We find the quantum critical 

points and their dependence on the anisotropy of the aforementioned field-induced quantum phase 
' transitions. The spin-fiop phase corresponds to the spontaneous breaking of Z2 symmetry. We use 

^ ' the numerical density matrix renormalization group and analytic spin wave theory to find the phase 

diagram of the model. The energy gap, sublattice magnetization, and total magnetization parallel 
' and perpendicular to the magnetic field are also calculated. The elementary excitation spectrums 

are obtained via the spin wave theory in the three different regimes depending on the strength of 
the magnetic field. 
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PACS numbers: 75.10.Jm, 75.50.Gg, 75.30.Ds, 64.70.Tg 



INTRODUCTION 



Quantum ferrimagnets are a general class of strongly correlated magnetism, which have attracted much interest 
' O . in experimental as well as theoretical investigations. Examples of such realizations are the bimetallic molecular 
magnets like CuMn(S2C202)2(H20)3-4.5H20 and numerous bimetallic chain compounds which have been synthesized 
systematicalljii^. In these materials, the unit cell of the magnetic system is composed of two spins, the smaller one 
\s a — 1/2 and the larger one (p) is changed from 1/2 to 5/2. The magnetic and thermodynamic properties of these 
models arc different from the homogeneous spin counterparts. For instance, the one dimensional mixed-spin model 
represents a ferromagnetic behavior for the low temperature regime while a crossover appears to the antiferromagnetic 
behavior as temperature increases^^— . The crossover can be explained in terms of the two elementary excitations 
. where the lower one has the ferromagnetic nature and a gapped spectrum above it with antiferromagnetic property^. 
^\ ' Moreover, the mixed spin models have shown interesting behavior for the quasi one dimensional lattices (ferrimagnetic 
ladders). Despite that the two-leg spin-1/2 ladder is gapful, representing a Haldane phase, the two-leg (mixed spin) 
ferrimagnet is always gapless with the ferromagnetic nature in the low energy spectrum. However, a special kind of 
, dimcrization can drive the ferrimagnetic ladder to a gapped phase^i^^. 
y—{ • The presence of a longitudinal magnetic field preserves the U(l) symmetry of the XXZ interactions and creates a 
T-H I nonzero magnetization plateau in a one-dimensional ferrimagnet for small magnetic fields in addition to the saturation 
• plateau for large magnetic fieldsii"— . The former plateau corresponds to the opening of the Zeeman energy gap 
which removes the high degeneracy of the ground state subspace. The ferrimagnets on ladder geometry present a rich 
structure of plateaus depending on the ratio and dimcrization of exchange couplingsi^. In both one-dimensional and 
two-leg ferrimagnets the magnetization plateaus can be understood in terms of the Oshikawa, Yamanaka and Affleck 
(OYA) argumen1>i^ because the longitudinal magnetic field commutes with the rest of the Hamiltonian and the models 
have U(l) symmetry. However, the situation is different when a transverse magnetic field is applied on the system, 
because the transverse field does not commute with the XXZ interaction and breaks the U(l) symmetry of the model. 
The onset of a transverse field develops an energy gap in a spin-1/2 chain which initiates an antiferromagnetic order 
perpendicular to the field directior>i^"— . The ordered phase is a spin-flop phase because of nonzero magnetization in 
the field direction; however, there is no magnetization plateau even in the gapped phased. The lack of U(l) symmetry 
prohibits the use of the OYA argument, thus prompts the question of a magnetization plateau and the presence of an 
energy gap2i in the spectrum. 

The structure of the paper is as follows. First we study the anisotropic ferrimagnetic chain in the presence of a 
transverse magnetic field by using the density matrix renormalization group (DMRG)^ and exact diagonalization 
Lanczos methods. The energy gap, sublattice magnetization, and total magnetization in both parallel and perpendic- 
ular to the field direction are presented in Sec. [Hi We further address the energy gap behavior versus the magnetic 
field and the magnetization plateau. The phase diagram of the model is also presented in the same section. We then 
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FIG. 1: The energy gap versus the transverse magnetic field. Difi'erent plots belong to various values of the anisotropy parameter 
A — 0.0, 0.5.1.0. Inset: The scaling of gap versus {2N)~^ for A = 0.0 at two critical points hd and hc2 confirms the vanishing 
of gap at these points while its scaling at ft = 2.0 verifies a finite gap in the thermodynamic limit (A'' ^ oo). 



use an analytical tool, the spin wave theory (SWT), to obtain the low energy excitation spectrum of the model in Sec. 
mil The SWT is applied in three different regions depending on the strength of the magnetic field. The qualitative 
behavior of the model is explained in terms of SWT and the magnetization is compared with DMRG results. The 
results of SWT help to explain the energy gap behavior of DMRG data. We finally summarize our results in Sec. IIVI 
where we put together both quantitative DMRG and qualitative SWT results to analyze the different phases of the 
model in the presence of a transverse magnetic field. 



II. DENSITY MATRIX RENORMALIZATION GROUP RESULTS 



We have implemented the numerical DMRG technique to study the magnetic properties of the anisotropic ferri- 
magnetic spin- (1/2, 1) chain in the presence of a transverse magnetic field given by the Hamiltonian ([T]): 

N 

H = J^[af pf + a\p\ + afpt^i + ^fP^Vi + ^{<P\ + - K^l + pf )], (1) 

1=1 

where erf (p.f ) represents the a-component of spin operators at site i for spin amplitude a = 1/2 (p = 1). The 
antiferromagnetic exchange coupling is J > 0, the anisotropy is defined by A, and h is proportional to the strength 
of the transverse magnetic field. 

The DMRG computations have been done on an open chain of length 108 spins (A^ = 54 unit cells) and the number 
of states kept in each step of DMRG is 300 < m < 500. We have also studied the chains with larger lengths (up to 
N = 100) and observed no significant changes on the data of magnetization and staggered magnetization within 5 
digits of accuracy. 

The energy gap is defined as the difference between the first excited state energy and the ground state energy. It 
shows whether the model is gapless or gapful depending on its zero or nonzero value, respectively. Using the DMRG 
computations, we have plotted in Fig. [T]the energy gap of the model versus the transverse magnetic field for different 
values of anisotropy parameter, A = 0, 0.5, 1.0. All plots show a gapped phase for small values of the magnetic field, 
h < hci{A), and a paramagnetic gapped phase for h > /ic2(A). The gap vanishes at two critical points, h = /ici(A) and 
h — ^c2(A). The isotropic case (A = 1) remains gapless in the intermediate region /id (A) < h < /ic2(A), while the 
anisotropic case (A 7^ 1) is gapful. The gap behaves differently for various A in the small- field and intermediate- field 
gapped phase. 

In the isotropic case A = 1 , the behavior of gap versus h can be explained in terms of the elementary excitations of 
the model. For A = 1, the U(l) symmetry of the model is restored and the magnetic field operator commutes with the 
rest of the Hamiltonian. Thus, the energy spectrum for A = 1 is expressed in terms of the spectrum at ft, = plus a 
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FIG. 2: (a) The a;-component sublattice magnetization versus a transverse magnetic field for both a = 1/2 and p = 1 spins 
and their sum as the unit cell magnetization in a::-direction for anisotropy parameter A = 0. (b) The i-component unit 
cell magnetization versus a transverse magnetic field for different anisotropics A — 0.0,0.5,1.0. The dashed line shows the 
saturation value at Mx = 3/2. 



shift of energy which depends on h. At h = the model has SU(2) symmetry and the ground state is a ferromagnetic 
state with total spin Sq = N{\p — a\), which is highly degenerate and the lowest ferromagnetic spectrum is a gapless 
one, namely v~{k) [see Eq. (fT7|) ]. An antiferromagnetic spectrum {v^ik)) exists above the ferromagnetic one, and 
the lowest state of the antiferromagnetic spectrum has total spin Saf = Sq + 1 with a finite gap 2J\p— (t|, measured 
from the ground state. Upon adding a commuting magnetic field to the ferrimagnctic chain the symmetry is lowered 
to U(l) and the energy levels arc affected by a Zeemann term, i.e., —hS^. For the magnetic fields smaller than h 
(which will be defined later), the Zeeman energy gain of the ground state is larger than all of the other states in 
the ferromagnetic spectrum; thus the ground state remains robust, and the first excited state is the first state in the 
ferromagnetic spectrum (with energy J{v~{k) + /i)), which leads to the energy gap equal to Jh. This explanation 
remains valid until the gain of the Zeeman term of the lowest state of the antiferromagnetic spectrum (J(i^"*"(fc) — h)) 
dominates the gain of the first excited state in the ferromagnetic spectrum. It defines h by the following equation: 

v~{k) + h = v+{k) -h, (2) 

which gives h = \p — a\ within linear approximation of SWT (from which both ^^{k) will be derived in the next 
sections). At this point, the first excited state is the lowest state of the antiferromagnetic spectrum. Thus, the energy 
gap behaves as {2J\p — a\ — Jh) before it vanishes a.t h = hci{A = 1) = 2\p — cr|. The linear increasing behavior 
for small fields and then linear decreasing of the energy gap are clear in the DMRG data for A = 1, shown in Fig. 
[TJ Although the DMRG values for h and hci have some discrepancies with those obtained by linear SWT, the SWT 
gives the qualitative behavior correctly. 

The energy gap of the anisotropic Hamiltonian (A 7^ 1) is defined as Ei — Eq for < h < hd and h > /ic2, where 
El is the first excited state energy and Eq is the ground state energy. However, the ground state becomes degenerate 
{El = Eq) for hci < h < hc2, where the energy gap is the difference between the second excited state energy and 
the ground state one, E2 — Eq. For small magnetic fields the scaling behavior of the energy gap can be explained 
using the quasi-particle excitations of the model as ft. — > 0. The leading term of quasi-particle excitations for very 
small magnetic fields {h — )- 0) gives the scaling of energy gap as Vh, for A 7^ 1 [in the weak field SWT, Eq.dlHl)]. In 
a similar manner, the leading term of the strong field SWT [Eq. ([2T|) ] leads to linear dependence of the gap on the 
magnetic field in the paramagneic phase which explains very well the behavior in Fig. [TJ The linear dependence of 
gap versus the magnetic field for h > hc2 is confirmed by the DMRG numerical data for any isotropics. 

We have plotted the energy gap versus (2iV)~^ in the inset of Fig. [T]to observe its finite size scaling (where 2N is 
the total number of spins). We have implemented both the Lanczos and DMRG algorithms to calculate the energy 
gap for A = 0. We have plotted the minimum value of gap which occurs at hd and hc2 versus (2A^)~^ which clearly 
shows that the gap vanishes in the thermodynamic limit (TV — ^ 00). It suggests that both hd and hc2 correspond to 
quantum critical points. The different magnetization characteristic confirms that a quantum phase transition occurs 
at both hd and hc2 (see Fig. [2]). Wc have also plotted the energy gap for h = 2.0 to justify that the gap of the 
intermediate region is finite in the thermodynamic limit. 

We have also plotted the x-componcnt magnetization of each sublattice in Fig. ^{a.) for ferrimagnctic spin-(l/2, 1) 
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FIG. 3: Unit cell magnetization (Mx) versus the anisotropy parameter (A) for some low magnetic field values (h). Our plots 
justify the plateau only for A = 1. The dashed line represents Mj; = 0.5 (the plateau value). 



chain with A = versus h employing the DMRG technique. The total magnetization has been plotted in Fig. [ll-(b) 
for different values of anisotropy, A = 0, 0.5, 1.0. To calculate the magnetization we have considered those spins which 
are far from the open ends of the chain to avoid the finite size boundary conditions. In this respect, ten spins have 
been neglected from each side of the open chain and the magnetization has been averaged over the rest of spins. 
Figure [2]-(b) shows the possibility of two plateaus in the magnetization along the field direction. For the isotropic 
case (A = 1), it can be explained in terms of the OYA argument^^. According to this argument, n{S — m) = integer, 
where n is the periodicity of the ground state, S the total spin of unit cell, and m a possible magnetization plateau 
of the unit cell, the one-dimensional spin-(l/2, 1) chain can show two plateaus at m = 1/2 and 3/2. However, for 
A 7^ 1 the axial symmetry of the model is broken by the transverse magnetic field, and the OYA argument is not 
applicable. Thus, more investigations is required to figure out the difference between the anisotropic (A 7^ 1) and 
isotropic (A = 1) cases. 

To get more knowledge on the behavior of magnetization for the anisotropic case, we have plotted the total mag- 
netization in the magnetic field direction (M^) versus the anisotropy parameter (A) for small magnetic field values 
in Fig. [3] The plots have been shown for those values of the magnetic field which seems to exhibit the magnetization 
plateaus. Figure [3] clearly verifies that the magnetization plateau only exists for the isotropic case, while there is no 
plateau for A 7^ 1. The magnetization per unit cell (Mx) in the direction of magnetic field {h) is given by 



1 dEo 
N~dh' 



(3) 



where Eq is the ground state energy. The above relation for a gapped phase simply states that if the ground 
state energy is linear in the magnetic field (Eq cx h), the magnetization will be constant, (the presence of plateau); 
otherwise the magnetization will depend on the magnetic field, (the absence of plateau). Let write the HaiTiiltonian 
'AS H = Hq — hHi where Hq is the XXZ interacting part and hHi is the magnetic field part. In the presence of 
U(l) symmetry (A = 1) the interacting and the magnetic field parts commute [Ho,Hi] = 0. Thus, Eq is a linear 
function of h which leads to the emergence of a magnetization plateau when the energy gap is nonzero. This agrees 
with the OYA statement. However, the transverse magnetic field breaks the U(l) symmetry in the anisotropic case 
(A 7^ 1) and [Hq^Hi] 7^ 0. Therefore, the ground state energy depends on h non-linearly which gives a change of 
magnetization when h varies, i.e. the lack of magnetization plateau even if a finite energy gap exists. 

Although the above general explanation is applied to the strong magnetic field regime the saturated plateau {M^ = 
1.5) can also be explained from another point of view. An eigenstate with full saturation is classified as a factorized 
state^^ in which all spins align in the direction of the magnetic field. As a general argument, it has been shown in Ref<^ 
that the full saturation for an anisotropic Heisenberg type interaction in the presence of a magnetic field takes place 
at a finite value of the magnetic field if the model is rotationally invariant around the field direction. Accordingly, 
the saturation at = 1.5 takes place only for the isotropic case (A = 1) and h > hc2- In the anisotropic case 
(A 7^ 1), the fully polarized plateau can take place for infinite strong magnetic field while the nearly saturated state, 
{Mx — 1.5), can be observed for large magnetic fields. To justify this argument we have plotted the x-component 
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FIG. 4: Schematic of spins' orientations in different phases of the anisotropic ferrimagnetic spin-(l/2, 1) chain in the presence 
of a transverse magnetic field. 



magnetization of each unit eell for different values of A in Fig. [2]-(b). It is clear that the magnetization in the field 
direction does not reach the saturation value of Mx = 1.5 for A = and 0.5, while it obviously touches its saturated 
value for A = 1 and /i > 3. 

The antiferromagnetic interactions between the spins in each unit cell make them to be antiparallel, which leads to 
the total a;-component magnetization = {ux + Px) — 0.5. This phase has been shown schematically in Fig. Ill-(l) 
where we have neglected the effects of small quantum fluctuations on the directions of the spins. The non-commuting 
transverse magnetic field opens a gap which is robust as long as h < hd- This (gapped) Neel phase corresponds 
to the first plateau at = 0.5 for A = 1 and a semi-plateau {M^ — 0.5) for A ^ 1. By further increasing h, 
the gap is closed at the first critical field /id (A) (for A = 0, hd ~ 1.6) where the magnetization starts to increase 
obviously. Further increasing of the magnetic field leads to a continuous change of the ground state property which 
gives a gradual change of the magnetization-Fig. |4]-(2-4). For strong magnetic field (/ic2(A = 0) > 2.4) the spins are 
nearly aligned in the direction of the magnetic field, the semi- plateau at ~ 1.5 [Fig. [2]-(a) and Fig. ID- (5)]. 

To get more insight on the ground state properties of the model, we have plotted the y-component spin expectation 
value versus the transverse magnetic field in Fig. [5]for A = 0. The magnetization in the y direction for both sublattice 
spins is zero for /i < 1.6 and h > 2.4; however, it becomes nonzero in the intermediate region 1.6 < h < 2.4. The values 
of the y component spins in the unit cell are equal, and their directions are opposite to each other, (ct^) = — (p^). It is 
surprising that for any value of the magnetic field 1.6 ^ h < 2.4 we get (cr^) = — (p^) whereas the spin magnitude on 
the sublattiees are different (0-7^ p). At the factorizing field, /i/ — 2.24 (which will be explained in the next section), 
where the condition crsinl^l = — psin|/3| should be satisfied, the mentioned relation is obtained (cr^) = — (p^). The 
staggered magnetization in the y direction, SAIy = (cr^ — p^), is nonzero for this region. Moreover, our numerical 
data verifies that the z component magnetization on both sublattiees is zero for any value of the magnetic field. 

Generally, let us consider the y component staggered magnetization as an order parameter, which is nonzero for 
/ici(A) < h < /ic2(A) and zero elsewhere. Nonzero SMy corresponds to a spontaneous breaking of Z2 symmetry. In 
fact, for any value of A the model has a Z2 symmetry which can be expressed by the parity operator P = (S)iafpf. This 
symmetry, which can also be considered as a tt rotation around the magnetic field direction (a;), leads to vanishing value 
for the z and y components of the spins. However, the symmetry is spontaneously broken for ft,ci(A) < h < /ic2(A) 
which selects one of the parity eigenkets to give nonzero sublattice magnetization in the y direction. 

III. SPIN WAVE ANALYSIS 

We have applied the spin wave theory to get more knowledge and a qualitative picture of the phase diagram. The 
SWT is a method to describe a spin model in terms of boson operators. The elementary excitations of the spin model 
are given by bosonic quasi-particles which are constructed on a given background. Based on this fact the SWT can be 
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FIG. 5: The y-component sublattice magnetization versus the transverse magnetic field for A = 0. This component is nonzero 
only in the intermediate phase 1.6 ^ h < 2.4. Moreover, the y-component spins are exactly equal and antiparallel for both 
sublattices. The staggered magnetization of the unit cell in the y direction is nonzero within the intermediate region. 



considered on different backgrounds to build up a bosonic system. Typically, a state in the Hamiltonian Hilbcrt space 
is considered as the background which is supposed to be the ground state within an approximation. However, there 
exists some spin models such as the isotropic antiferromagnetic Heisenberg spin-1/2 chain that do not have an ordered 
ground state, and thus the SWT fails to explain the properties of the model correctly. Therefore, the existence of an 
exact ground state is a good starting point to initiate a spin wave analysis. 

The ferrimagnetic chain both in the absence and presence of longitudinal magnetic fields has been studied by the 
SWT— iSiii^ii^. Although the Nccl state is not the exact ground state for a ferrimagnct in the presence of a longitudinal 
field, the SWT gives a good description of the model which justifies that the quantum fluctuations arc not strong 
enough to ruin the whole picture. It would be more interesting to initiate a spin wave theory based on an exact 
ground state for a ferrimagnct in the presence of a transverse magnetic field. According to Ref.— , the exact ground 
state of a general class of ferrimagnets can be found at the factorizing magnetic field, h = hf. This ground state is 
a factorized state, which is a perfect background to implement SWT. It gives a reliable analysis around h = hf (see 
next subsection). We will also study the SWT for small and large magnitudes of the magnetic field. Our analysis 
is limited to the linear spin wave approximation to get the magnetic properties of the spin-((T, p) ferrimagnets in the 
presence of a transverse magnetic field. 



A. SWT at h = hs 

Let us briefiy introduce the exact factorized ground state of a ferrimagnct in the presence of a magnetic field^^. 
The factorized ground state for the Hamiltonian of Eq. ([T|) can be written in the following form: 

i^o)= (g) Wi)w;). (4) 

where \a'^ and are the eigenstates of Gi ■ n[ and pj ■ h'- with the largest eigenvalues, respectively, with and 
being unit vectors pointing in polar angles (6*, </? = 0) and (/3,a = 0). A„ and Bp represent the two sublattices which 
contain the two different spins. The factorized state is called a bi-angle state, defined by the two angles {0,/3) and 
represents the ground state of the model at h = hf, where 

p + Aa 

cosp 



(T + Ap 

Vp2+CT2 + 2A^' 



(5) 



7 



and 

hf = 2 + 0-2 + 2Aap, 
ef = -{a^ + p^ + Aap), (6) 

with ey being the ground state energy per site at the factorizing field. 

To perform the spin wave analysis around h = hf, we first implement a rotation on the original Hamiltonian (H). 
The rotated Hamiltonian (H) is the result of rotations on all lattice points of H, and is given by the following relations 

H = D^'HD, 

D = (g) Dt{Q,e,0)DP[0,p,Q). (7) 

The rotation operator 

^"(0, /3, 0) = D{a = 0, 7 = 0) = D,{a)Dy{l3)D,{-f), 

is defined in terms of Euler angles, and a similar expression is considered for C(0, 6*, 0). 

In the rotated basis defined by {x\y',z') and {x" ,y" , z"), the bi-angle state becomes the fully polarized ground 
state of H. In the next step, the rotated Hamiltonian is bosonized via a Holstein-Primakoff (HP) transformation. 



aZ = \/2cr - a,|aj a^, erf = a - aja. 



2p-b]b,b,, pf^p-b]h,, (8) 

where a,i{a\) and bjib^-) are two types of annihilation (creation) boson operators, satisfying the commutation relations: 
[ai,a]] = 6ij, = Si^, [ai,6]] = and [a,;,6j] = 0. 

The Hamiltonian in the momentum (fc) space and in the linear spin wave approximation is written as 

El = NJapcos{/3 ~ 9) + NJh{pcos (3 + a cosO), 



Hi = jY,l^2^Acos'^{akbl + bkal) 



H2=J 



2tT cr 

]V(/i; - h) 



Vpsin/3 (feo + 



+^/a•sin^^(ao + aj) 

(9) 

where is the total number of spins in each sublattice. The unitary transformation that diagonalizes Hi is given by 

Xk = ak cos T]k - bksimjk, 

il^k = bk cos r]k + flfc s'mijk, (10) 

where Xk and ipk are the quasi-particle boson operators that preserve the bosonic commutation relations. In this 
representation, we obtain 

Hi^Y.'^u;-{k)xlxk+u;+{k)4i;k), (11) 
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in which are the quasi-particle excitation modes. The dispersion relations are given by 

+ 2A^/CTptan(277fc) cos | 



w*(fc) = JD+±J- 



v/1 + tan2(277fc) 



2AVpCT cos 



tan(27?fc) = ^ (12) 



in which 



- f(^±^) + Mjcos^±^cos/3) 

- A(g±p)+ '^ (cos/gzfccosg). (13) 

A shift on the zero momentum component of boson operators, defined by two constants t^, diagonahzes the full 
Hamiltonian; i.e., Xo — ^ Xo + ^~ 7 4'o ^ i'o + , where 

^/2N{h — hf) (-yo- sin 770 sinfl + cos 770 sin/3) 



2lo+ 

'2N{h — hf) (-y/CTCOS77o sinO — sin r/o sin /3) 



2uj- 



(14) 



The diagonalized Hamiltonian is given by 



k 

Egs = Ei+uj-f.-^ +uj+t^ 



+V2N J{hf — h) ^(v^ sin 6* cos 770 — ^/psin-qQ sin/3)t 



+ ('\/CTsin0sin7^o + y^cos 7/0 sin /3)t^ 

Ef = 2NJef (15) 

where Egs is the ground state energy which reduces to Ef at the factorizing field (/i/) (i.e., the energy of the exact 
bi-angle state). 

The magnetic properties of model ([T]) can be studied through the linear spin wave theory-Eq. ^T5\i. In Fig. |6]-(a), 
we have plotted the sublattice magnetization of the anisotropic ferrimagnetic spin- (1/2, 1) chain for A = 0.5. The x 
and y components of sublattice magnetization are nonzero [Fig. [B]-(a)]; however, the 2;-component of the sublattice 
magnetization is zero, denoting that the spins are located in the xy plane. It should be noted that the values of 
sublattice magnetization is exact at the factorizing field while it is approximately correct for the magnetic field close 
to the factorizing field. The a-component of total magnetization per unit cell is Ma = (p" -t-cr") and the corresponding 
staggered magnetization is defined SMa = (p" — cr"). In Fig. ini-(b), we have plotted the x and y components of total 
magnetization and staggered magnetization. Around the factorizing field the model has a considerable magnetization 
in the x direction and a staggered magnetization in the y direction, which identifies a spin-flop phase around the 
factorizing flcld. The model has a dual character i.e. it behaves like a ferromagnet in the x direction and like an 
antiferromagnet in the y direction, it is the result of two branches of excitations, Eq. (jl5|) . which are the origin of the 
existence of two dynamics in the model^iiS. We will discuss later the effects of magnetic field on the configuration of 
both spins in more details. 



B. SWT at weak and strong magnetic fields 

(a) Weak Field SWT 

In the SWT it is assumed that the ground state defines a particular classical direction for the spins. In the weak 
magnetic fields close to h = 0, we expect to have a Neel-ordered configuration. Therefore, we use the following 
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FIG. 6: (a) The sublattice magnetization, (b) Magnetization and staggered magnetization per unit cell of the anisotropic 
ferrimagnetic spin-(l/2, 1) chain versus transverse field, for A = 0.5. The factorized ground state is chosen as the background 
in the linear SWT. 



Holstein-PrimakofF (HP) transformations: 



'^t = al Y 2a' - alflj, cr,f = -cr + aja^, 

p+ = ^2p-b]b, b, , p^^p- b]b, . (16) 

In the linear spin wave approximation and within Fourier space representation, one can diagonalize the Hamiltonian 
which is given by 

H = Eo + J2{'^-ik)vlvk + iy+ik)wiwk}, (17) 

where 

Eo = -A^J(2o-p + p + a)- A^J/i(p-a) + i^(zy-(fc) + i.+ (fc)), 

fc 



,y±{k) = J\/2(p2 + s2 _ 2Apc7COs2 ^ ± Di), 
Di = y(p2-s2)2-4[A(p2 + s2)-ps(l+A2)]p(TCOs2| 

h h 
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FIG. 7: The sublattice magnetization, total magnetization and staggered magnetization per unit cell of the anisotropic ferri- 
magnetic spin-(l/2, 1) chain versus transverse field, for A = 0.5, when the Neel order is chosen as the background in the linear 
SWT. 



and Vj., w^(vfe, Wk) are bosonic quasi-particle creation (annihilation) operators. The procedure of the diagonahzation^S 
dictates that the bosonic Hamiltonian should be positive definite. This constraint implies that for |A| < 1 the 
amount of magnetic field obeys the condition h < 2{p — a), and for 1 < |A| < the magnetic field should be 

\h- p + a\< v^(p + (t)2 -4/9crA2. 

Let us consider the special case oi {a = = 1) and A = 0.5. The Hamiltonian of this system (in the linear SWT 
approximation) is positive definite only for magnetic fields smaller than ho — 1. Accepting this condition we have 
plotted in Fig. [7]the sublattices field-induced magnetization, the total magnetization and the staggered magnetization 
per cell of the whole chain versus transverse field h. It shows that for < /i < 1 the model is affected slightly by the 
transverse magnetic field. In other words, the staggered magnetization in the x direction is close to its maximum value 
(the Nccl ordered state). The quantum fluctuations for < h < 1 are not strong enough to change the magnetization 
from its zero field value. However, upon reaching ho ~ 1 the quantum fluctuations are suddenly increased so that 
they destroy the ordered state completely. Thus within this linear SWT, the first critical field is hf^'^ = = 1 
and for an arbitrary (cr, p)-ferrimagnetic chain it becomes hf]^^ = 2\p — tj|. Although hf^"^ does not depend on 
the anisotropy parameter A and is slightly different from the DMRG results (Fig. [2]), the linear SWT describes the 
elementary excitations of the model well. 

(b) Strong Field SWT 

For the strong magnetic fields the ground state is ordered in the direction of the magnetic field. The fully polarized 
ground state in which all spins are aligned in the field direction is used as the background for initiating the SWT. In 
this case the following HP transformation is implemented for the spin operators 



CTj = cr — aiOi, 

p^j^p-b]b,. (19) 
The diagonalizcd Hamiltonian in terms of the Fourier space representation and within the linear SWT is 

H = Eo + J2{^'ik)V^Vk + n+{k)wlWk}, (20) 

k 
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FIG. 8: The magnetization of sublattices, the total magnetization, and the staggered magnetization per unit cell of an anisotropic 
ferrimagnetic spin-(l/2, 1) chain versus transverse field and for A = 0.5 and when the background in the linear SWT is the 
field-induced fully polarized state. 



where 



Eo = 7VJ(2ap + p + cr)- Arj/i(p + cr + l) + i^(f)-(fc)+l)+(fc)), 



D2 

P 



(21) 



and ,Wl{Vk,Wk) are bosonic quasi-particle creation (annihilation) operators. The condition to have a positive 
definite bosonic Hamiltonian implies that for |A| < 1 the amotmt of the magnetic field shotild be larger than 2(p + a) 
and for |A| > 1 the magnetic field should be larger than p + a + y/{p — cr)^ + ApaA'^. 

Again we consider the special case of (cr = = 1) and A = 0.5. The Hamiltonian of this system in the linear 
SWT approximation is positive definite only for a magnetic field larger than hf^"^ = 3. The magnetization of each 
sublattice, the total field-induced magnetization, and the staggered magnetization per unit cell are plotted in Fig. 
[51 For ft, > 3, the model is in the polarized phase. We have already shown in Refi ^'^i^^ that the full saturation only 
happens for the isotropic case A = 1. Thus the model possesses an upper critical field hc2 = 3 for A = 1. The 
comparison with DMRG results shows that hf2^'^ = 3 is the true value, which is the consequence of weak quantum 
fluctuations for the strong field regimes. For A 7^ 1, the fully saturated state appears at infinite magnetic field. It 
can be understood simply by imposing = = /3 in Eq. ([5]) which can be fulfilled only for A = 1 in the Hamiltonian 
given by Eq. ([I}. In general, the full saturation occurs at a finite magnetic field if the model has the U(l) symmetry 
around the direction of the magnetic field. 

Let us discuss qualitatively the effects of a non commuting transverse magnetic field on the phase diagram of the 
anisotropic ferrimagnetic spin-(l/2, 1) chain. The SWT gives two branches of quasi-particle excitations for each of 
the small, intermediate and large magnetic field regions. At zero magnetic field the lower branch is gapless with 
ferromagnetic nature while the upper one is gapped with antiferromagnetic signature. A nonzero magnetic field opens 
a gap in the ferromagnetic branch which remains robust for h < hci- Moreover, the staggered magnetization in the 
field direction is close to its maximum value which implies a Neel phase. At /i = hci a quantum phase transition 
from the Neel phase to the spin-flop phase takes place where the staggered magnetization perpendicular to the held 
direction becomes nonzero. The quasi-particle excitations for the spin-flop phase are given by uj^(k). In the spin-flop 
phase {hci < h < hc2) an entanglement phase transition occurs at h ~ hf where the quantum correlations become 
independent for h < hf and h > hf. The increase of magnetic field causes the second quantum phase transition at 
h = hc2 to a nearly polarized state in the field direction. The excitations in the field induced polarized phase {h > hc2) 
are gapful given by 17^ (fc), where the gap is proportional to the magnetic field. 
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IV. SUMMARY AND DISCUSSION 

The ground state phase diagram of the anisotropic ferrimagnetic (ct, p) chain in the presence of a non commuting 
transverse magnetic field has been studied. The general picture has been obtained within the spin wave approximation. 
We have applied three schemes of linear spin wave approximation to find the magnetic phase diagram of the anisotropic 
ferrimagnetic spin- (cr, p) chain with anisotropy parameter A and in the presence of the transverse magnetic field (h). 
The spin wave approximation has been applied close to /i = (weak fields), h = hf (intermediate regime), and 
h ^ hf (strong fields), where hf is the factorizing magnetic field. The ground state is known exactly aX h = hf 
as a product of single spin states. We have studied the magnetization in the field direction. There is a plateau 
at Mx = 0.5 for isotropic case where the ground state energy is linear in magnetic field while no plateau observed 
for the anisotropic cases. However, the magnetization along the magnetic field changes slightly as long as ft. < hci 
and its value is Mx — 0.5, which motivates to recognize it as a Neel phase . The model exhibits a quantum phase 
transition at /i = hd from the Nccl phase to (i) a spin-flop phase for A 7^ 1, (ii) a gapless Luttinger liquid for 
A = li^iH. The magnetization evolves in the spin-flop phase when the magnetic field is increased. The spin-flop phase 
contains the factorizing field {h = hf) where an entanglement phase transition takes place and quantum correlations 
vanish. Further increase of the magnetic field leads to a polarized phase which resembles a plateau at the saturated 
magnetization in the field direction. However, it will be fully saturated only for A = 1 (the presence of a rotational 
symmetry around the magnetic field) which is represented by a quantum phase transition at a finite value /ic2. The 
validity domain of spin wave analysis were introduced and it was shown that the corresponding results were in good 
agreement with the DMRG numerical computations. 

To get more accurate values on the magnetization process of spin-(l/2, 1) fcrrimagnct, we have also plotted in Fig. [5] 
the DMRG data of the x- and y-componcnt staggered magnetization in addition to the x-component magnetization of 
unit cell versus the transverse magnetic field for A = 0. The magnetization curve has been divided to five regions which 
has been labeled in fig. 21 fig. [HI and also in Table. H] Region-(l) is defined by the Neel phase for < /i < hd — 1.6 
where both Mx and SMx are nearly constant while SMy is zero. The spin-fiop (gapped) phase, hd < h < hc2, 
where a nonzero SMy sets up can be distinguished to three parts, namely regions- (2-4). For hd <h< 1.9 which is 
labeled region-(2) we observe (cr^) < and (p"^) > 0. It is a spin-fiop phase which is called spin-fiop (I) in Table. HI 
Rcgion-(3) is defined at /i ~ 1.9 where the projection of smaller spin along the magnetic field becomes zero, (cr^) — 0, 
i.e. Mx = SMx- The rest, 1.9 < /i < /ic2 — 2.4, labeled by rcgion-(4) where (cr^) > and {p'') > is called spin-flop 
(H). The region-(5) is the polarized phase along the direction of magnetic field, i.e. Mx — 1.5 and SMy = 0. It is 
observed from Fig. ^{a) that the component of smaller spin in the direction of the magnetic field is affected strongly 
by the magnetic field while the corresponding component for the larger one is almost constant. 

The spin-fiop (I) is a characteristic behavior of XXZ ferrimagnets in the presence of transverse magnetic field 
because the spin component of the smaller spin along the magnetic field is opposite to the field direction ((cr^) < 0) 
while the spin-flop (II) is similar to the corresponding phase of the homogeneous XXZ spin chain in the presence of 
transverse magnetic flcld {{(J^) > 0)^i^. In the anisotropic ferrimagnetic chain the transverse field first develops a 
Neel phase and a field-induced quantum phase transition leads to a spin-flop phase. Moreover, the Z2 symmetry is 
spontaneously broken for small-field region in the homogeneous spin chain while it will be broken in the intermediate 
fields /ici(A) < h < /ic2(A) for ferrimagnets. A summary of different properties of the homogenous XXZ spin 1/2 
chain and the corresponding (1/2, 1) ferrimagnet both for isotropic and anisotropic cases is presented in Table. |TT1 



TABLE I: Different configurations of the ground state of the ferrimagnetic spin-(l/2, 1) chain with A = in the presence of a 
transverse magnetic field. 



Region 


h 


Phase 


Order parameters 


(1) 


< /i< 1.6 


Neel 


AU = 1/2, SMy = 


(2) 


l.G<h< 1.9 


Spin-Flop(I) 


((T^-> < 0,SMy > 


(3) 


h ~ 1.9 


Spin-Flop 


{a^) = 0,SAIy > 


(4) 


1.9 < /i< 2.4 


Spin-Flop(II) 


(<T^-> > 0,SMy > 


(5) 


h > 2.4 


Nearly Polarized 


Mx 3/2, SMy = 



It is also interesting to mention that the low energy effective Hamiltonian of the anisotropic spin- (1/2, 1) chain in 
the presence of a transverse magnetic field can be represented by the fully anisotropic (XYZ) spin-1/2 Heisenberg 
chain in an applied field (though we do not report such calculations in this paper) . This helps to get more knowledge 
from the results on the effective model^l. However, both spin wave approximation and DMRG results show that the 
model has two nearly constant magnetization in the presence of transverse magnetic field, the small-field plateau at 
Mx — 0.5 for h < /ici(A) and the saturated Mx — 1.5 for large fields {h > ftc2(A)). The general behavior is the 
same for any value of the anisotropy parameter (A); however, the critical fields /ici(A) and /ic2(A) depend on A. For 
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h (transverse magnetic field) 



FIG. 9: The a;-component magnetization, x- and y-components staggered magnetization versus the transverse field for a 
ferrimagnetic spin-(l/2, 1) chain. Effects of the magnetic field on the spins of each sublattice are divided into five different 
regions. 



instance, /ici(A = 0.5) ~ 1.8 and /ic2(A = 0.5) ~ 2.6. 

TABLE II: Different ground state phases are classified for the heterogeneous spin-(l/2, 1) XXZ ferrimagnet along with the 
homogeneous spin 1/2 XXZ antiferromagnet. The comparision between isotropic and anisotropic cases in the presence of the 
transverse magnetic field [h) is presented. The magnetization per unit cell is m. The ferrimagnet has two critical points hd 
and hc2 while the homogeneous antiferromagnet has a critical point at h^. 



Spin 


Region 


Isotropic case (A = 1) 




Anisotropic case (A j^- 1) 


(1/2,1) 


< h < hci 


Gapped Neel, plateau at m = 1/2 




Gapped Neel, no plateau 


(1/2,1) 


hci < h < hc2 


Gapless Luttinger liquid, no plateau 




Gapped spin-flop, no plateau 


(1/2,1) 


h > hc2 


Gapped paramagnet, plateau at m = 


3/2 


Gapped paramagnet, no plateau 


1/2 
1/2 


0<h<hc 
h > he 


Gapless spin-fluid, no plateau 
Gapped paramagnet, plateau at m = 


1/2 


Gapped spin-flop, no plateau 
Gapped paramagnet, no plateau 



The magnetization process can also be viewed as a non-unitary evolution of the system. The entanglement of a 
pure state (ground state in our ease) is conserved under local unitary operational^. For the ferrimagnetic spin- (1/2, 1) 
chain, the entanglement of the system is decreased by increasing the magnetic field for h < hf. The entanglement 
vanishes at h = hf where the ground state is given by a tensor product state. This is an entanglement phase transition. 
It is thus concluded that the effect of magnetic field is a non-unitary evolution of the ground state. 
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